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We investigate the heat flow in the parametric quantum 
pump. Using the time dependent scattering matrix theory, we 
have developed a general theory for the pumped heat current 
at finite pumping amplitude and frequency. We have applied 
our theory to a double barrier structure and studied pumped 
heat current in both the weak and strong pumping regimes as 
different system parameters vary. By comparing the pumped 
heat current and the power of Joule heat generated in the sys- 
tem, we found that the double barrier structure can function 
as an optimal pump in the strong pumping regime. 
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The physics of adiabatifi-nuantuni pump has attracted 
great attention recentlyErta. The quantum pump is real- 
ized by varying the geometric parameters of the quantum 
dot, by which the current is generated. At the same time, 
heat current is also produced and accompanied with the 
dissipation. Avron et alij have given a lower bound for 
the dissipation in a quantum channel which is defined as 
the difference between the heat current and the power 
of Joule heat. The pump is optinaal if the heat current 
equals to the power of Joule heatEj. As a result, the op- 
timal pump is noiseless and charge transported is-quan- 
tized. In a recent paper, Moskalets and ButtikeinJ have 
also considered the dissipation in an adiabatic quantum 
pump. The heat current and the noise has been formu- 
lated in terms of a parametric emissivity matrix. The 
theory of Moskalets and Buttiker is in the weak pump- 
ing regime (quadratic order in pumping amplitude) and 
can go to the finite frequency if one goes beyond the hy- 
pothesis of instant scattering. In this paper, we develop 
a general theory for the heat current which is valid for 
finite pumping amplitude and finite frequency. This al- 
lows us to study the heat current in both weak and strong 
pumping regimes. Our theory is Ipai^fd on the time de- 
pendent scattering matrix theoryn^pEl and goes beyond 
the instant scattering hypothesis. We have applied our 
theory to a double barrier structure and studied pumped 
heat current as different system parameters such as Fermi 
energy, pumping amplitude, and phase difference vary. 
In the weak pumping regime, the pumped heat current 
increases quadratically as pumping amplitude increases. 
The dependence becomes linear in the strong pumping 
regime. In the strong pumping regime, the pumped heat 
current shows strong nonlineality as a function of phase 
difference between two pumping potentials. The heat 
current for single pump has also been studied. We found 
that the amplitude of the pumped heat current for the 
single pump is of the same order as that of two pump- 



ing potentials. Recently there is a concern of whether 
a genuine optimaLpump with nonvanishing transmission 
coefficient existsEil. In this paper, we give a nontrivial 
example of optimal pump. By comparing the pumped 
heat current and the power of Joule heat generated in 
the system, we found that the double barrier structure is 
an optimal pump in the strong pumping regime. 

We start with the general defeiition for the heat cur- 
rent in scattering matrix theorycJ, 
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Iq,a = lim -r- / dt < > 

At->oo At Jq 
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where the heat current operator is Iq^a — lE,a — Eple^a 
and < ... > denotes the quantum average. Here lE,a is 
the energy current operator given by 

= ~t[dtblit)bo,{t) - dtaiit)a^it)] (2) 

and Ie,a is the electric current operator, 

ie,c. ^ bUt)bait) - aUt)aa{t) (3) 

where the operators ba and Ua are annihilation opera- 
tors for the outgoing and incoming carriers in the lead a. 
They are related by the scattering matrix. 



(4) 



where the time-dependence of the scattering matrix is 
due to the slowly time- varying pumping potential X{t). 
The distribution fitmction can be obtained by taking the 
quantum aver age, c3 
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where aa{E) is the Fourier transform of aa{t) and f{E) 
is the Fermi distribution function. For the purpose of 
presentation, we calculate the energy current first and 
the electric current can be calculated in a similar fashion. 
From Eqs.(|), (|), and the energy current is given 

by 
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dt I dtidt2^^Sai3{t,ti) 
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where f{t) = J {dE/2-K) exp{-iEt)f{E). Now we will fo- 
cus on the first term (denoted as /jj^l^) in Eq.(@). After 
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changing of the variable to = {ti + t2)/2 and t = ti — ^2 
and using thp, foUowing Wigner transform for the scat- 
tering matrix 



s{t,t') 



Eq.(0) becomes, 
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[dts'^pi^2, -) + «e2S„/3(e2, — 7)] 
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Changing the variables again to ri — t ^ to and t' = 
(t + to)/2 and integrating over ri, we obtain 

/At 
Jt'drdeY.sMe,t' + r/mr) 

[(l/2)dt>sl^{e, t' ~ t/4) + zesl^ic t' - t/4)] e^^ (9) 

To get the heat current up to w^, it is enough to expand 
Saf3 up to the second order in r. We obtain 
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dt / dE[2dEf - Edlf] 
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where the zeroth order term in r in Eq. (|^) has been can- 
celed by the second term in Eq.(^ and Tp is the period 
of the pumping cycle. Including the electric current, the 
heat current is given by 
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Y,dtS^p{E,t)dtslp{E,t) 



dt / dE[2dEf -{E- EF)dlf] 
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dt / dE{E - EF)dEf 



2TrTp Jq 
Y,Mdtsi^iE,t).s„p{E,t)] 
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Note that the second term in Eq.(|ll[) vanishos at zero 
temperature. In the adiabatic regime, we havea dtSap = 
J^il^XiSapdtXi + dj^.SapdtXi + ...] where X = dX/dt. 



Up to the order up' , we can neglect the contribution from 
dx.Saf}- At zero temperature, Eq.(pj]) becomes. 



^Xj^^ap 9tXi dfXj 



(12) 



It is straightforward to obtain the heat current at higher 
order in frequency. To do that, we have to keep the higher 
order term in the expansion of r in Eq.(01 and include 
the contribution of higher order derivativesau such as X. 
For instance, expanding Eq.(||) to the third order in r 
and including the contribution from the electric current, 
we obtain (up to uj'^), 



li^l = I ' dt I dEdlfY Im 



(13) 



There are also third order corrections from Eq. (|ll|) where 
the term d-^ Sap should be kept. We now consider the 
limiting case of Eq.(^2|) when the pumping amplitude is 
small. For two probe pumping: Xi{t) — sin(a;t) and 
X2{t) = X2 sin(a;t + (f), we see that the lowest order 
pumping amplitude in Eq. ( p^ is quadratic by neglecting 
X dependence in dxSgg. In this case, it is straightfor- 
ward to show that Eq. (O) is reduced to 
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2cos(l) XiX2j2Mdx,Sg0dx,sl^)] (14) 



which agreea23 with the result of Ref. 

We now apply our formula Eq.(|12|) to a one- 
dimensional quantum structure which modeled by double 
barrier potential U{x) = XiS(x + a) + X26{x — a) where 
2a is the well width. For this system Ihe Green's func- 
tion G{x,x') can be calculated exactlycj. With G{x,x') 
we can ^calculate scattering matrix from the Fisher-Lee 
relationBj Sap = —Sap + ihvG{xa, xp), with v the elec- 
tron velocity in the lead. The adiabatic pump that we 
consider is operated by changing barrier heights adia- 
batically and periodically: Xi = Vio + Vp sin(a;t) and 
X2 — V2o + Vp sm{ujt+(l)). This can be achieved by micro- 
fabricating metallic gates at the barrier region and apply- 
ing a time dependent gate potential. Since the pumped 
current is proportional to o;^, we will set uj = l-4or con- 
venience. Finally the unit is set by ?i = 2m = lES. 

We first study the pumped heat current with two 
pumping potentials. Fig.l depicts the pumped heat cur- 
rent exiting from left lead versus Fermi energy for the 
symmetric barriers at small pumping amplitude. We 
have also plotted the transmission coefficient (solid line) 
versus Fermi energy for comparison. The physical pi|e-| 
ture of heat flow suggested by Moskalets and ButtikerE2l 
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is the following: as electron scattered by the oscillating 
barriers or scatterers, the absorption of energy quantum 
huj creates an electron-hole pair. The flow of electron- 
hole pair leads to the heat transfer. In the symmetric 
case, we have Iq^L = Iq,R- We see that the heat current 
is peaked at the resonant levels and is clearly propor- 
tional to the density of states of the scattering region. 
At phase difference (p = 3tt/4 (short-dashed line), the 
lineshape of Iq is approximately Lorentzian similar to 
that of transmission coefhcient. As (j> decreases, the line- 
shape is broadened and deviates from Lorentzian consid- 
erably. For large pumping amplitude (see inset of Fig. 1), 
we see that there is significant heat current even in the 
off resonant case (when Fermi enegy is not in line with 
the resonant level in the static case when pumping po- 
tential is off). This is because in the strong pumping 
regime, the instantaneous resonant level oscillates with a 
large amplitude and hence can generate heat current in 
a broad range of energy. Up to the order of w^, the heat 
currents Iq^L and Iq^R are all positive and flow from the 
scatterers to the reservoir. It is no longer true when high 
order frequency contribution is included. Fig. 2 shows the 
heat current as a function of phase difference for different 
pumping amplitudes. In general, there are two extreme 
points for the heat current at = and tt, where 
the former corresponds to maximum in heat current and 
latter for the minimum. In the weak pumping regime, 
the sinusoidal behavior is seen. In the strong pumping 
regime, however, we see signiflcant derivation from the 
sinusoidal behavior. As a result, when increasing from 
= to TT, the heat current drops not as fast as in the 
weak pumping regime away from (p = tt. However, near 
(/> = TT, the heat current decreases much faster. Fig. 3 
displays the heat current as a function of pumping am- 
plitude. The general behaviors of heat current at differ- 
ent (j) are similar. We see that initially, the heat current 
increases quadratically with pumping amplitude in the 
weak pumping regime and then quickly approaches to the 
linear regime in the strong pumping regime. This behav- 
ior persists for single pumping potential. In the inset, we 
present the result for asymmetric barriers. Here ?;io = 4, 
V20 — 6, and Vp — 0.5, Iq^L is about twice as large as Iq^R- 
This is reasonable because it is easier for heat current to 
tunnel through the lower barrier. In general, the electric 
pumped current is linear in frequency for two pumping 
potentials whereas for single pumping potential, it must 
be zero up to the first order in frequencyaila. The heat 
current is different. For both single pump and two poten- 
tial pump, the heat currents are of order w^. In the weak 
pumping limit, there is a simjple relationship between the 
following three heat currentsE2l: 1^/2 — II + Ir- where 
is the heat current of two pumping potentials with 
phase difference (j) — vr/2, is the heat current of sin- 
gle pump due to the left oscillating barrier, and Ij^ heat 
current due to the right barrier. In the strong pumping 
regime, the scattering matrix depends on both Xi and 
X2 in a nonlinear fashion, this simple relation is no longer 
valid. Denoting the ratio t = /7r/2/(-^L + Ir)^ we found 



that at small pumping amplitude, r ~ 1. As one in- 
creases the pumping amplitude the ratio r decreases and 
quickly saturates around r ~ 0.8. This is understand- 
able since the behavior of 1,^/2: II, and Ir are similar 
as one increases the pumping amplitude, one obtains the 
nearly constant ratio in the linear pumping amplitude 
dependence regime (see Fig. 3). 

Finally, we study the dissipation of the pump by com- 
paring the heat current and the power of Joule heat pro- 
duced by the electric currentcj. We found that the heat 
current is always larger than the corresponding power 
of Joule heat as predicted by Ref. In Fig. 4, we 

plot the pumped charge per cycle and the ratio between 
the power of Joule heat and the heat current as a func- 
tion of pumping amplitude Vp. We see that the pumped 
charge increases as Vp increases. In the large Vp limit, the 
pumpe d . ch arge will eventually reach the maximum value 
Q = ellil'tZl. We see from Fig. 4 that the ratio also ap- 
proaches to one just like the pumped charge. Therefore, 
we conclude that the double barrier pump we studied can 
be optimal in very strong pumping limit. 

In summary, we have developed a general theory for 
the pumped heat current at finite pumping amplitude 
and frequency which allows us to investigate the pumped 
heat current in both weak and strong pumping regimes 
as different system parameters vary. As the pumping am- 
plitude varies, we observed the crossover of pumped heat 
current from quadratic dependence in the weak pump- 
ing regime to linear dependence in the strong pumping 
regime. In the strong pumping regime, the pumped heat 
current shows strong nonlineality as a function of phase 
difference between two pumping potentials. For the sin- 
gle pump, we found that amplitude of pumped heat cur- 
rent is of the same order as that of two pumping poten- 
tials. Finally, our numerical results show that the double 
barrier structure we examined can be an optimal pump 
in the strong pumping regime. 
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FIG. 1. The transmission coefficient (solid line) and heat 
current as a function of Fermi energy at diff'erent phase differ- 
ences between two pumping potentials: = (dotted line), 
4> = tt/2 (dot-dashed line), (ji = 37r/4 (short-dashed line). The 
heat current for a single pump (long-dashed line). For illus- 
trating purpose, the heat current has been multiplied by 10/3. 
Other parameters: uio = 1120 = 80 and Vp = 4. Inset: the heat 
current as a function of Fermi energy at different pumping 
amplitude: Up = 5 (dotted line), Vp = 10 (dot-dashed line), 
Vp = 20 (dashed line). The transmission coefficient is also 
shown (solid line). The heat current has been normalized to 
one. The scaling factors are: 8/3 for Vp = 5, 4/3 for Vp = 10, 
and 2/3 for Vp — 20. Other parameters: wio = ^20 ~ 80, 
(f> = 0. 



FIG. 2. The heat current as a function of phase differences 
(I)/tt at different pumping amplitudes: Vp — 1 (solid line), 
Vp — 20 (dotted line), Vp = 40 (dotted-dash line). The heat 
current for a single pump (long-dashed line). For illustrating 
purpose, the heat current has been multiplied by a factor 
of 200, 20/3, 10/3, respectively for Vp = 1,20,40. Other 



parameters: iiio 



■ "20 



80 and Ep = 2.4069. 



FIG. 3. The heat current as a function of pumping ampli- 
tude at different phase difference: <j) = (solid line), (p = Tr/2 
(dotted line), (p — 3tt/4: (dotted-dash line). The heat current 
for a single pump (dashed line). The heat current has been 
multiplied by a factor of 10/3. Other parameters are the same 
as the inset of Fig. 2. Inset: Iq^L (solid line) and Iq^R (dot- 
ted line) versus pumping amplitude for asymmetric barriers. 
Here (p = and Ef = 1.65. The scaling factor is 1/3. 



FIG. 4. The pumped charge per cycle (solid line), the 
ratio between the power of Joule heat and the heat current 
(dotted line) as a function of pumping amplitude. The system 
parameters are t;io ~ V20 = 80, — S/pi/i, and Ef = 2.4069. 
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